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SUMMARY 

A fully second-order closure model for turbulent reacting flows Is sug- 
gested based on Favre statistics. For diffusion flames the local thermodynamic 
state Is related to a single conserved scalar. The properties of pressure 
fluctuations are analyzed for turbulent flows with fluctuating density. Clos- 
ure models for pressure correlations are discussed and modeled transport 
equations for Reynolds stresses, turbulent kinetic energy dissipation, density- 
velocity correlations, scalar moments and dissipation are presented and solved, 
together with the mean equations for momentum and mixture fraction. Solutions 
of these equations are compared with the experimental data for high heat 
release free mixing layers of fluorine and hydrogen In a nitrogen diluent. 


INTRODUCTION 

A test case for the evaluation of a fully second order closure model of 
turbulent reacting flows with moderate to large heat release has been provided 
by the recent experiments of Hermanson et al. (ref. 1) and Mungal and Dimotakis 
(ref. 2). The dynamic field consists of a two-dimensional mixing layer with 
gas phase free streams carrying hydrogen In a nitrogen diluent In one stream 
and fluorine In a nitrogen diluent In the other. The reaction Is 

H 2 + F 2 *-2 HF 

and Is highly exothermic. The experimental results with the highest heat 
release correspond to fluorine concentration of up to 6 percent and hydrogen 
concentration of up to 24 percent, with a maximum adiabatic flame temperature 
rise of 940 K. 

An attempt will be made to assess the possible Improvement In the pre- 
diction of the above flow field by accounting for both velocity and scalar 
fluctuation time scales and by Inclusion of the chemical heat release effects 
In the modeling of the Reynold stresses, scalar fluxes, and dissipation 
equations. 


Closure Model 

The closure model consists of three parts: 

(1) The thermo-chemical model based on the shifting equilibrium flame 
model (ref. 3). This model requires a fast rate of chemical reaction and 
relates the value of all thermodynanlc variables to a conserved scalar. There 
are several scalar variables which are conserved In a chemical reaction and can 
be used as a basis for describing the mixing In a nonpremlxed reacting flow. 

The mass fraction of a given element Is such a variable. A normalized 





conserved scalar Is referred to as the mixture fraction. Hence the density 
p, temperature T, and composition are local functions of the mixture 
fraction f. 

(2) The turbulence or moment model consisting of a set of closed Favre 
averaged equations of order one and two. Including the equations for all 
Reynolds stress components, scalar fluxes, dissipation rates of turbulence 
kinetic energy, and the mixture fraction variance. 

(3) The coupling model which relates the thermo-chemical model and the 
turbulence model through the probability density function, p.d.f., of the mix- 
ture fraction. In general, the form of p.d.f. will depend on the flow condi- 
tions and will be coupled with the chemical heat release. If constraints are 
Imposed by solving the transport equations of the mean and the variance of the 
mixture fraction, the errors that can be made using an arbitrary form of the 
p.d.f. are limited. A Favre Beta function will be used here to represent the 
p.d.f. of the mixture fraction. The means and higher moments of the thermo- 
dynamic variables at any point may be obtained from their local relationship 
with the mixture fraction obtained from the thermo-chemical model and the 
p.d.f. of the mixture fraction. 


NOMENCLATURE 

D velocity divergence 

f mixture fraction f = f + f" 

fs Stoichiometric value of the mixture fraction 

g^ gravitational acceleration vector 

k turbulent kinetic energy 

p pressure, p = p" + p 1 

£ Favre probability density function 

R ratio of turbulence velocity to scalar time scales 

T absolute temperature 

T ^ adiabatic flame temperature 

T ^ ambient temperature 

U longitudinal velocity component 

~ w 

velocity vector, + u^ 

V transverse velocity component 

x position vector 

mass fraction of the k^ species In the mixture 
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'0.5 


r 



p 


T 


<P 


transverse location of the point at which (ft - - U,,) = 0.5 

dlffuslvlty (same for all species) 

Integral thickness 

vortlclty or maximum slope thickness 

shear layer 1 percent thickness of the mean temperature profile 

turbulent kinetic energy dissipation rate 

scalar variance dissipation rate 

density of the mixture 

time scale 

equivalence ratio 


( ),( )' conventional averaging 

(~),( ) " density weighted averaging 

( )^.( ) 2 high speed and low speed sides 

TURBULENCE MODEL 

Second-order closure methods are believed to be the optimum level of 
closure and have been successfully applied to constant density flows. Includ- 
ing buoyancy effects (refs. 4 to 6). The presence of the density fluctuations 
In the variable density or combusting flows results In far more complex equa- 
tions for the second moments. Application of density weighted (Favre) averag- 
ing to variable density flows results In a set of equations which are similar 
to those of constant density flows. This extends the use of the well developed 
constant density second-order closure models to the variable density case and 
makes available well tested model expressions for many terms In the moment 
equations (ref. 7). This Is, however, not sufficient for the complete closure 
of the density weighted moment equations, since new correlations Involving 
density and divergence of velocity appear In the set. The exact equations show 
that the variation of density Is felt via three mechanisms: (1) the variation 

of the mean density and the correlation of density with velocity; (2) the non- 
zero divergence of velocity; and (3) the correlations involving pressure which 
are related to density via mass, energy, and state equations. These correla- 
tions appear In the transport equations of the Reynolds stresses and scalar 
fluxes and play an Important physical role. 


Pressure Equation 

The Instantaneous pressure equation follows from mass and momentum balance 
as 



(pU^Uj) + g^p + a 
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(1) 


3 



where denotes the Newtonian stress tensor. In contrast to the con- 

stant density case, the type of this equation for known velocity Is hyperbolic. 
It describes the change of pressure due to hydromechanical motion and acoustic 
propagation of waves. The density term however contains not only acoustic 
variations of pressure but contributions which are due to mixing and chemical 
reaction. Expanding the first and second terms of above equation and utilizing 
the Instantaneous mass and momentum equations, the above equation can be 
rewritten as, 


ap . -,. lWl * - yy - oy - - yy tj * yy 


where the substantial derivative Is 


D t - a t * V, 


( 2 ) 


and 


D , 3,11, - - ; y 

The last two terms In equation (2) represent the molecular viscous effects on 
pressure and are Ignored In comparison to other terms In high Reynolds number 
flows. The second term In equation (2) results In an Integral equation of the 
Fredholem type of second-order. It can be shown that the Iterated kernals for 
this equation exist, thus a solution can be obtained. This solution Is rather 
complicated, but It produces an Integral form with a kernal modified by the 
density fluctuations. This term will have negligible affect at high Reynolds 
numbers due to the lack of correlation. The first term on the right-hand side 
of equation (2) Is the well known hydrodynamic source of pressure change, see 
Chou (ref. 8). This would be the only term appearing In the case of constant 
density flows. In the context of the thermochemical model for diffusion flames 
Introduced above, density Is a local function of the mixture fraction. It Is 
however advantageous to consider p-1 , specific volume, Instead of p. 

Therefore 

D = p df. p D^f 

The pressure equation can be recast in terms of mixture fraction, 

Ap a -pd^UjdjU^ - pD^ ^p d f p ^fj (3) 

In pure mixing of two gases of different molecular weights the variation of the 
specific volume as a function of the mixture fraction Is close to linear, with 
Its second derivative equal to zero. However In diffusion flames this relation 
Is strongly nonlinear near the stoichiometric value of the mixture fraction. 

To the left and right of this value p-T ( f ) Is nearly linear with different 
slopes. If this local relation Is simplified to partially linear with a jump 
In the slope at the stoichiometric value of the mixture fraction, then the 
terms Involving the first derivative of the specific volume reflect the effect 
of mixing, with 
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(4) 


d f/3 _1 . a(l - H (f - f s )) 

and the terms Involving the second derivative of the specific volume reflect 
the effect of chemical heat release, with 

d = -a6(f - f $ ) (5) 

where H(f - f s ) denotes the Heaviside function, 6(f - f s ) the Dirac function 
as generalized derivative of H(f) and f s is the stoichiometric value of 
the mixture fraction. Utilizing the specific volume, the equation of conserva- 
tion of mixture fraction and considering the limiting case of high Reynolds 
number and Kolmogorov's hypothesis of local Isotropy (refs. 9 and 10) the 
pressure equation for diffusion flames is then reduced to 

Ap » -p3 1 Uj3 J ajU 1 + 2 P 2 E 2 d 2 f (ln p) d 2 f p _1 (6) 

where E = ra^fd^f. Based on the above model and equation (5) the second 
term on the right hand side of equation (6) is recognized as the heat release 
contribution to the pressure equation. This term behaves like a Dirac delta 
function and its magnitude is determined by the location and movement of the 
Instantaneous flame front, f(x,t) = f s , since only near the flame front it 
has a large contribution and is small everywhere else. This simplified pres- 
sure equation analysis indicates that the fluctuations of the pressure contain 
the contributions of mixing and reaction in addition to velocity fluctuations. 
Janicka and Lumley (ref. 11) have argued that the mixing part can be neglected 
at high Reynolds numbers in pure mixing. An order of magnitude analysis shows 
that in the case of chemical heat release the mixing contribution can also be 
neglected. 

Solution of the Poisson equation, (eq. (6)), for pressure can be obtained 
by application of Green's theorem. Thus far from walls, where the surface 
Integral may be neglected 


P(x) - h 00 + P2(x) 

where 

P l (x) = U J j"J pa j U i a i U j |x - x' | dx l dx 2 dx 3 (7) 

and 

P 2 (X) = " 2^ J J J p ( ra i f3 1 f ) d ff (ln d ff p 1 |x - 1 x. ' | dx l dx 2 dx 3 

(8) 


Second-Order Closure 

The complete set of equations constituting a closed system of Favre- 
averaged moments including all density variation effects will be discussed 
briefly. 
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The mean velocity, at high Reynolds number, satisfies 


3 4 pU i u jj = -v - pg i 

No closure assumption Is required at this level of closure. 
The modeled Reynolds stress equation Is given as 


(9) 


3 1 pU l U 1 U j = 


— ii it ^ — ii ii ^ | 2 — 

pU 1 U J 3 1 U 1 + pU ] U^ 3 ] Uj I - - peS^ 


+ 3 


1 



U,U 3 
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l ( ii _ ii _ 

+ z jp' u i a j P + p' u j a ^ p 




( 10 ) 


The Kolmogorov assumption of local Isotropy at high Reynolds numbers Is used 
here, so that the dissipation tensor Is Isotropic with Its trace equal to twice 
the turbulence kinetic energy dissipation, e. The diffusion model of Daly and 
Harlow (ref. 12) Is chosen over the more complicated model suggested by 
Hanjallc and Launder (ref. 13) due to their equivalent performances In constant 
density flows. Consideration of the exact solution of the pressure equation 
given by equations (7) and (8), suggests the following decomposition of the 
velocity-pressure gradient correlation 
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The first term on the right-hand side of the above equation Is denoted by 
q}j and represents the correlation of the velocity fluctuations with the 
hydrodynamics part of the pressure fluctuation gradient. This would be the 
only term present for the Incompressible flow case and contains the pure 
turbulence Interaction "return to Isotropy," and mean strain rate Interaction 
"fast response" contributions. The closure model of Launder, Reece, and Rodl 
(ref. 4) Is used for this term. In density weighted form this term Is given as 


Q 1j = ~ C 1 p k \Vj " 


3 4 1j k 


C 2 + 8 
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"ij " 3 6 1j W 
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u^u^S^Uj + UjU^a^U^J and ir = -u^Uj3^U^ 


I II II ^ II II ^ \ II II 

^ = - (u^u^U-i + UjU^a^U^ I and k = 1/2 u^u^. 
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The term represents the heat release part of the pressure- 

strain rate correlation. Insight Into modeling of this term can be gained by 
careful analysis of the pressure equation's solution (eq. (8)). Considering 
just the heat release portion of the pressure solution we have 



( 12 ) 


Due to the Dirac delta function behavior of the second derivative of the 
specific volume near the stoichiometric value of the mixture fraction, this 
correlation Is determined by the location and movement of the Instantaneous 
flame front, f(x,t) = f s , l.e., only near the flame front does the second 
derivative of the specific volume contribute to the pressure equation. There- 
fore, the Integral representation of the heat release contribution Indicates 
dependence on the probability of the flame sheet being present at a given 
location. On the other hand, the probability of mixture fraction being In an 
Interval about the stoichiometric value, £(f s ;x) df, Is zero for pure oxidizer, 
where f = 0, and grows quickly with fuel concentration. It reaches Its maxi- 
mum at the flame front region and declines to zero In the pure fuel region, 
where f = 1. This suggests that f s J2i) Is large and positive on the 

entrainment side of the flame and becomes small and negative past the flame 
front. Considering the dimensional properties of equation (12) It becomes 
clear that there are two time scales affecting this term. One Is the scalar 
fluctuation time scale defined by the scalar variance and Its dissipation, the 
other one Is the dynamic or velocity fluctuation time scale, which Is related 
to the turbulence large eddy structures. If tensorlal properties and con- 
sistency with the limits of zero fluctuations and pure mixing without reaction 
are required, the following possible closure model emerges: 



(13) 


where R = 2 ef/f" 2 k/e denotes the ratio of turbulence time scale to scaler 

fluctuations time scale and cf = ra^Pa^f" Is the rate of dissipation 
of the scalar variance. Beguler, Dekeyser, and Launder (ref. 14) showed from 
a survey of several nonlsothermal , nonreacting turbulent shear flows that the 
time scale ratio was approximately equal to 2.0. The Implied constancy of R 
Is used In most work aimed at calculating scalar variances. However based on 
the evidence to the contrary provided by Warhaft and Lumley (ref. 15) an exact 
transport equation for the scalar dissipation Is derived and modeled, so that 
R can be determined at each point. 

The constant Cri can be estimated as 

C R1 ~ C * 8* d ff ln p 


where C* Is a measure of correlation among the terms In the Integral 
(eq. (12)) and should have a value much less than one due to the degree of 
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statistical Independence Implied In the modeling of this term, d^f In p 
Is a chemically dependent term that varies from zero (for nonreacting cases) 
to the order of 100 for reactions with large heat release. Therefore the value 
of C R i depends on the energy release In the flow field. For the hydrogen 
- air flame of Dibble, Kollmann, and Schefer (ref. 16) a value of 1.25 Is 
recommended (ref. 10). 


The rate of turbulence kinetic energy dissipation determines the dynamic 
length and time scales and satisfies the following modeled equation 



u.u.a. c 
1 J 1 i 


- c e i p 


k u i Vj°i - C c2' 
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(14) 


where the last term contains the effect of chemical heat release on dissipation 
This complex process Is tentatively taken Into account by analogy with the 
corresponding term (eq. (13)) for the normal stresses. 


Density - Velocity correlations appear both directly and through the 
modeled terms In the Reynolds stress transport equation (10). Derivation of 
exact equations and the closure of the equations are discussed In detail by 
Farshchl (ref. 10). The modeled transport equations for these correlations Is 
given by 
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— II II II 


II II / 11 


+ 3j^C s tpUjU 1 3 1 p'u 1 j- 2C s TU j U 1 3 1 ^p , U i J 3 jp + ~ 3^P 


H I ii 


^V /2 n 2 ^' 2 


" C RU0 Vj ^ + C 1RU i «*' U 1 - C RF k f R p »1 


+ c 2 ruPp' U J (15) 

t Is the turbulent time scale, which Is traditionally set equal to the 
turbulence time scale k/e. However It can also be the time scal e of scalar 
fluctuations or a combination of the two time scales (ref. 17). p' 2 Is the 

variance of density and Is obtained directly using the thermo-chemical model 
and the appropriate weighted Integration of the mixture fraction p.d.f. 

The transport equation of the mixture fraction's mean value requires no 
modeling and Is given as 
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(16) 


a^pi^T) 




where the molecular diffusion effects are neglected In comparison with the 
turbulence diffusion at high Reynolds numbers. 


Next a transport equation for the scalar flux vector, u f " , Is derived 
and modeled (ref. 6). The scalar pressure gradient terms (pressure scrambling) 
are modeled using the same line of argument used for velocity-pressure gradient 
correlations In the Reynolds stress equation. The modeled form of this equa- 
tion Is given as 





II It ^ 

* oVjy 


I — II II II 

+ a, c Tpu,u,a,u.f" 
j V c ^ j 1 1 1 , 


+ = pf'^P - C lc 1 pujf" + C 2 c pUjf“ 3 j U 1 


( n o\l/2 p -if » — ^ 

f £ ) R^ ~ £ ~~ — ^(p£(f s ;x)) (17) 

7 p 

p'f" appearing In the above equation can be obtained from the p.d.f. of the 
mixture fraction. 

Mixture fraction variance appears In the Beta function representing the 
probability density of the mixture fraction. For high Reynolds number flows 
the modeled mixture fraction variance equation Is expressed as 
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1 



— II ^ 

-2pU^f"3^f + 3^ 


f II II II 

. C c2 TU l“jV 




(18) 


The variance of the mixture fraction Is dissipated by molecular action at small 
scales. This dissipation rate Is also the key to determination of scalar 
fluctuation time scale. The exact form of the scalar dissipation equation and 
Its modeling Is discussed In details by Farshchl et al . (ref. 18). However, 
regrouping some of the exact terms and applying an order of magnitude analysis 

-1 2 

and Ignoring terms of order Re. and less, where R = k /vc Is the turbu- 

L G t 

lence Reynolds number, Farshchl (ref. 10) shows thatfof high heat release cases 
there Is a source along with a sink of the scalar dissipation. The chemical 
source term for diffusion flames Is given by 


d ff'( r W’) 2 * c p3 „E( f s ;x) .f i 

f £ 

with c p 3 = 1.0. The completely modeled scalar dissipation equation Is then 
given by: 
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(19) 


All constants used In above equations are summarized in table I, 


APPLICATION AND COMPARISON WITH MEASUREMENTS 

Hermanson's (ref. 19) experimental studies of the effect of heat release 
in a planar, gaseous reacting mixing layer formed between free streams contain- 
ing hydrogen and fluorine in nitrogen diluents are used for comparison. The 
density ratio of the free streams is kept equal to one by using as diluent a 
mixture of nitrogen and a small amount of helium, on the fluorine side, and a 
mixture of nitrogen with small amount of argon on the hydrogen side. The 
nominal high speed, hydrogen containing side, flow velocity Is 22 m/s and a 
free-stream speed ratio of Up/U] = 0.4 is maintained. The high speed 
free-stream turbulence level was measured to be about 2/3 percent. The measur- 
ing station was positioned 0.457 m downstream of the splitter plate trailing 
edge. The Reynolds number at the measuring station based on the high speed 
free-stream velocity, the downstream distance and the cold free-stream kine- 
matic viscosity was Re x = 6xl0 5 . The Reynolds number based on high speed 
boundary layer momentum thickness just upstream of the trailing edge was 
estimated to be about 240. 

Two cases of moderate and high heat release are considered, corresponding 
to two equivalence ratios (ref. 2): 

(1) «p = 1, corresponding to fluorine concentration of 6 percent and 
hydrogen concentration of 6 percent, with a maximum adiabatic flame temperature 
of 860 K. 

(2) (p = 1/4, corresponding to fluorine concentration of 6 percent and 
hydrogen concentration of 24 percent, with a maximum adiabatic flame tempera- 
ture of 1240 K. Thermodynamic tables relating density and temperature to the 
mixture fraction are obtained using CEC83 (ref. 20) computer program. 


Numerical Method 

The numerical method used for this calculation is a modified version of 
the Patankar and Spalding method (ref. 21) where all the first order moments 
are defined on node points, and ail higher order moments are defined at 
midpoints. The given Initial velocity profiles have Blasius form and agree 
with the boundary layer momentum thickness value reportedly the experiment. 
The Initial dissipation rate, e, and the shear stress, u"v", are determined 
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utilizing the turbulence equilibrium assumption, a simple gradient flux model 
with molecular kinematic viscosity, and the assigned velocity profile. Next 
the well known turbulence viscosity model, 


„ t - 0.09p jj- 

with laminar viscosity Is used to determine Initial turbulent kinetic energy 
profile. The normal stresses are equl-partitloned . The Initial mean mixture 
fraction has a step profile and all other turbulent correlations are set to 
zero. The longitudinal mean pressure gradient must be prescribed analytically 
In terms of the gradient of the outer boundaries velocities and Is zero here. 
The transverse pressure gradient cannot be prescribed and Is approximated by a 
simplified version of the mean transverse momentum equation as 


9 y P = -3 y (pv'' 2 ) 

The upper bound of the marching step of the main procedure Is determined at 
each forward position with respect to the value of a characteristic thickness 
of the mixing layer. To solve a typical problem with this 13 equation turbu- 
lence model over a distance of 0.9 m with 80 node point In the cross direction, 
the time needed on an IBM 370 Is about 30 min CPU. The p.d.f. Integration for 
the calculation of mean density, temperature, and other mixture fraction 
moments at each node takes 2/3 of above time. 


Discussion of the Results ‘ 

The Initial prediction of the homogeneous cold flow mixing layer resulted 
In a mixing layer thickness 35 percent below Its experimental value at the 
axial measuring station. To validate the model constants and the method used 
here and due to the lack of velocity correlations measurements by Hermanson 
(ref. 19) It was decided to use Browand and Latlgo's (ref. 22) experimental 
study. Two measures of thickness of the turbulent mixing region are employed. 
First Is the Integral thickness defined as 


S t u i - U(y) ][U(y) - U ] dy 
1 ( AU; -oo i 


Second Is the vortlclty thickness defined as 


6 w = 


JML 


^ 3 y U ^max 


Figure 1 shows that there Is about 16 percent difference between the experi- 
mental and calculated values of the. normalized Integral thickness. The axial 
rate of change of vortlclty thickness Is experimentally estimated to be between 
0.15 and 0.17 X, where X = (U-| - U 2 )/(U-| + U 2 ). The calculated value, d x 6 w = 
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0.148 X, Is at most 13 percent less than the experimental value. Comparison 
of normalized mean velocity and turbulence shear stress are presented In fig- 
ures 2, 3, and 4. They show good agreement with experiment except near the 
high speed side where the calculation shows a faster approach to the outer 
values . 

Physically the spread of the mixing layer downstream of the splitter plate 
Is determined by two different mechanisms. The spreading rate of the mixing 
layer at dynamical equilibrium stage Is one mechanism. This Is reasonably well 
predicted by the calculation. The other and more crucial mechanism, which Is 
not predicted by the calculation, Involves the transition from attached bound- 
ary layers on the splitter plate to free mixing layer developing downstream of 
It. If the boundary layers are laminar then there Is also a transition from 
laminar to turbulent flow In the developing region of the mixing layer. Birch 
(ref. 23) has attempted to model the developing region of the mixing layer by 
a multi-length scale gradient flux turbulence model. However, present calcula- 
tions Indicate that a single-length scale second-order model along with sub- 
stantially reduced value of the Initial turbulent kinetic energy dissipation 
would result In an Increase of the spreading In the Initial developing region. 
Even though the final spreading rate In the similarity region Is the same as 
before, the thickness of the mixing layer Is larger by up to 25 percent. The 
major short coming of this treatment Is that the amount of Initial reduction 
of e Is not the same for all flow conditions and Is a function of velocity 
ratio and character of the boundary layers on the splitter plate. Since In the 
study of heat release effects only relative values of the spreading are of 
Importance the basic model without any modification can be used. 

To examine the contribution of heat release terms Introduced by pressure 
fluctuation correlations and those appearing In the scalar dissipation equation 
two sets of model constants are examined: 

(I) Constants Cr-j , Cr 2 . Crp, Cq 3 , and C p 3 are all set to zero 

(II) The values given In table I are used for above constants 
The rest of the constants are kept the same In both cases. 

The comparison of mean axial velocities (fig. 5), and turbulent shear 
stresses (fig. 6), obtained for the large heat release case, <p = 1/4, utiliz- 
ing models (I) and (II) Indicate that there Is a little difference between the 
two models and the Improvement obtained by the application of model (II) Is 
almost negligible. This means that the heat release portion of the pressure 
fluctuation correlations has little effect on the turbulent velocity field. 

This result can also be shown by the fact that at each level of heat release, 
models (I) and (II) predict very close values for the spreading rate of the 
mixing layer. The spreading rate, however, decreases with the heat release. 
Spreading rate Is defined by d ( yo T - yo.9)/ dx » ^- e -* the rate of spread 
of the distance between the points at which U equals 0 . 1 ( U-j - U 2 ) + U 2 
and 0.9(U-| - U 2 ) + U2. Table II summarizes the spreading rates at different 
levels of heat release. 

The moderate heat release case causes about 10 percent reduction In the 
spreading rate, however the Increase In the heat release has only produced 
about 3 percent further reduction. The spreading rate levels off as the heat 
release Is Increased. This leveling corresponds with the leveling of the mean 
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density reduction reported by Hermanson (ref. 19) and leads to the conclusion 
that probably the single most Important cause of the spreading rate reduction 
with heat release Is the reduction of density and the subsequent reduction of 
the turbulent shear stress, -pu"v". Turbulent shear stresses normalized by 
free-stream density and velocity difference are compared In figure 7. There 
Is a considerable reduction of shear stress between the cold and moderate heat 
release cases. Further Increase of the heat release, however, results In a 
small reduction In the value of shear stress. 

Hermanson (ref. 19) reported a monotonic decrease of the vortlclty thick- 
ness, 6 W , with heat release without any Influence of the equivalence ratio 
on this trend. The present equilibrium flame formulation of the chemical 
reaction Is however dependent on the stoichiometric value of the mixture frac- 
tion and therefore on the equivalence ratio, <p. For 6 percent fluorine and 
<p = 1 , the stoichiometric value of the mixture fraction Is about 0.5 and for 
<p = 1/4 It Is about 0.2. It Is In the neighborhood of this value that the 
mixture density reaches Its minimum value. Figures 8 and 9 show the differ- 
ences In the density variance distributions and the mean density minimum value 
location for the two cases of tp = 1 and 1/4. The movement of the mean den- 
sity minimum value causes the shift In the profile of turbulent shear stress 
and the location of Its peak value as displayed In figure 7. Since the veloc- 
ity profile Is dependent on the gradient of the turbulent shear stress and not 
Its absolute magnitude, the vortlclty thickness predicted for the case of mod- 
erate heat release, tp = 1 , Is smaller than the value predicted for the high 
heat release case, <p = 1/4. Hermanson's (ref. 14) experimental results for 
6 percent fluorene and tp = 1 also show some scatter In vortlclty thickness 
measurements with some values that are smaller than those measured for the 
higher heat release case, «p = 1/4. These values are compared with calculated 
results In table III. 

Considering the general trend of the experimental vortlclty thickness 
behavior and assuming that above scatter Is within acceptable measurement 
accuracy bounds, one concludes that the above discrepancy between the predicted 
and measured results and the dependence on the stoichiometric value of the 
mixture fraction can be resolved by accounting for Intermlttency effects by 
such methods used by Kent and Bllger (ref. 24). Byggstoyl and Kollmann 
(ref. 25) have presented a model for the prediction of Intermlttency and 
Initial exploratory calculations with such models point to better predictions. 

Figure 10 compares the predictions of model (I) and (II) for the ratio of 
scalar to velocity field time scales, Indicating that the Influence of modifi- 
cations Included In model (II) Is larger on the scalar field than the velocity 
field. Even though the difference In density predictions (fig. 11), Is not 
large, the density variance and the mixture fraction variance predictions 
(figs. 12 and 13), show as much as 36 percent jump In their values with model 
(II) near the stoichiometric value of the mixture fraction. This change In the 
magnitude of the mixture fraction variance causes more than 10 percent reduc- 
tion In the peak value of the predicted temperature which corresponds very well 
with the experimental measurement, as shown In figure 14. The transverse 
coordinate In this figure Is normalized by shear layer 1 percent thickness, 
which Is defined as the transverse width of the layer at which the mean tem- 
perature Is 1 percent of the maximum mean temperature rise. The shift between 
the predicted and measured temperature profiles Is symmetric on both sides and 
Is expected to be a consequence of a fixed horizontal high speed upper side 
wall and diverging low speed lower side wall In the experiment. 
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Comparison with the limited experimental results available Indicate that 
the modeling approach used here Is a step In the right direction. However, 
comprehensive measurements of the scalar field quantities, such as mixture 
fraction and density variances and scalar-velocity correlations, are needed to 
examine and fully evaluate the contribution of the model presented here. 


CONCLUSIONS 

A fully second-order closure model for reacting turbulent flows have been 
extended to Include the effects of heat release on he turbulence flow field. 
This was done by analyzing the exact pressure equation of a turbulent nonpre- 
mlxed flame. The scalar field time scale has also been directly calculated by 
developing an equation for the dissipation of scalar variance which also 
Includes the heat release effects. Two models, (I) not Including heat release 
terms and (II) Including the heat release terms were compared and Indicate 
that: (a) The velocity field predictions are very similar In both models, 
correctly predicting the reduction In mixing layer spread rate with heat 
release. The scalar field results of model (II) are In better agreement with 
the experimental data yielding a 10 percent Improvement In temperature pre- 
dictions. (b) The time scale ratio Is sensitive to density fluctuations. This 
points to a correct Inclusion of heat release effects In the equation of scalar 
variance dissipation. Nevertheless, It would be premature to draw any conclu- 
sion about universality of such a closure and the values of the constants used 
here. Further experiments are needed to support or to Improve present 
assumptions and modeling. 
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TABLE 

I. - 

MODEL 

CONSTANTS 

c s - 

c i 

C 2 

C R1 

c 

G 

C el C c2 

0.2 

1.6 

0.43 

1 .25 

0.15 

1.45 1.90 

r 

c 

c 

c 

C 

C 

c3 

R2 

1 RU 

2RU 

RF 

RUD 

1 .0 

1 .25 

4.0 

0.5 

1.25 

1.0 

C 

c 

C 1C 

C 2C 

C C2 

C D1 

C D2 C D3 

0.18 

°pl 

1.8 

3.2 

C P2 

0.8 

0.5 

C P3 

1.0 

1.18 

2.2 

0.8 1.0 


TABLE II. - SPREADING RATES OF THE MIXING 
LAYER AT DIFFERENT LEVELS OF 
HEAT RELEASE 


Heat release level 

d(yo.i - yo.9) /dx 

Cold flow 

0.0562 

Moderate heat release 

.0505 

(«p = 1 » T adb = 860 K) 


High heat release 

.0488 

(«p = 1.4, T adb = 1240 K) 



TABLE III. - COMPARISON OF NORMALIZED VORTICITY 

THICKNESS WITH MEASUREMENTS 

w wo 


[6 W0 Is the cold flow vortlclty thickness.] 



Experimental (ref. 14) 
measurements 

Present 

calculation 

<p « 1 moderate 

0.7645 

0.8814 

heat release 

.8066 

.8814 


.9067 

.8814 

<p = 1/4 high 

.8826 

.9627 

heat release 
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Figure 1. - Normalized local integral thickness as a function 
of downstream distance for X = 0.695 and 6 tj = 0.46 mm. 



Figure 2. - Longitudinal mean velocity. 
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Figure 5. - Comparison of velocity prediction of models 
I and II with the experiment for the high heat release 
case, <p = 1/4. 



Figure 6. - Comparison of turbulent shear stress predic- 
tion of model I and II for high heat release case, <p - 1/4. 



COLD FLOW 



Figure 7. - Comparison of turbulent shear stress for 
different heat release values, using model I. 



Figure 8. - Comparison of density variances for different 
equivalence ratios, using model I. 
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